Dynamical Signatures of Edge-State Magnetism on Graphene Nanoribbons 
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We investigate the edge-state magnetism of graphene nanoribbons using projective quantum 
Monte Carlo simulations and a self-consistent mean-field approximation of the Hubbard model. 
The static magnetic correlations are found to be short ranged. Nevertheless, the correlation length 
increases with the width of the ribbon such that already for ribbons of moderate widths we observe a 
strong trend towards mean-field-type ferromagnetic correlations at a zigzag edge. These correlations 
are accompanied by a dominant low-energy peak in the local spectral function and we propose that 
this can be used to detect edge-state magnetism by scanning tunneling microscopy. The dynamic 
spin structure factor at the edge of a ribbon exhibits an approximately linearly dispersing collective 
magnonlike mode at low energies that decays into Stoner modes beyond the energy scale where it 
merges into the particle-hole continuum. 

PACS numbers: 71.10.Fd; 73.22.Pr; 75.40.Mg 
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Graphene is regarded as a promising candidate for fu- 
ture electronics. The electronic properties of nanodevices 
based on this material are in a large part governed by 
their edge structure. Besides their potential for nano- 
electronics and spintronics [H-Ill, graphene nanoribbons 
exhibit a remarkable interplay of low dimensionality, a bi- 
partite lattice and electron-electron interactions. One of 
the most fundamental predictions for graphene nanorib- 
bons (with zigzag edges) is the possibility of spontaneous 
edge-state magnetism, whereas bulk graphene is nonmag- 
netic Although electron-electron interactions are 
crucial for the emergence of edge-state magnetism, essen- 
tially all these predictions are based on mean-field type 
approximations. While such mean-field approximations 
were found - by comparing to numerical exact results - 
reliable at least for certain quantities and in the weak- 
coupling regime 13| on small honeycomb clusters, Lieb's 
theorem predicts nonet magnetization for a bipartite sys- 
tem at half filling 141] and the Mermin- Wagner theorem 
[ist actually forbids true long-range order for a ribbon 
with a fixed finite width. This is consistent with a sigma- 
model treatment which predicts a spin gap, i.e., a finite 
spin-spin correlation length for an even number of zigzag 
lines jl6j |. Furthermore, a numerical study indeed found 
a quite short correlation length in the case of two zigzag 
lines 
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In view of this puzzling situation, and given the grow- 
ing experimental effort toward probing the magnetism 
of graphene nanoribbons, we believe that there is need 
for an accurate numerical treatment of the Hubbard 
model description. Here we present projective ground- 



state quantum Monte Carlo simulations [3| and com par e 
them with a self-consistent mean- field theory (MFT) [l3| . 

For narrow ribbons we find indeed a finite correla- 
tion length; i.e., the static magnetism is an artifact of 
the mean-field approximation. However, with increas- 
ing width we observe a fast convergence of the quali- 
tative behavior of the correlation functions towards the 
MFT result corresponding to strong ferromagnetic or- 
dering tendencies. Since the ferromagnetic order is not 
static but subject to fluctuations, we next explore the 
dynamical properties of the nanoribbons, focusing both 
on local and momentum-resolved single-particle and spin 
spectral functions. We find that the single-particle spec- 
tral function on the edge exhibits a low-energy peak as 
a dynamic signature of the edge-state magnetism. The 
spin spectral functions reveal the onset of magnonlike ex- 
citations on top of the fluctuating magnetic background 
along the zigzag edge. We exhibit a linear contribution 
to their low-energy dispersion [?, 1^ 2^ and, even more 
interestinglyjthe damping of this mode into the Stoner 
continuum [20|. 

The ribbon geometry considered here is depicted in the 
inset of Fig. [T] We employ periodic boundary conditions 
at the armchair edge and open boundary conditions at 
the zigzag edge. L and W denote the length and the 
width of the ribbon, respectively. L is measured in units 
of lattice vector ai = (\/3, 0) and W is the number of 
zigzag legs of the ribbon. The carbon-carbon bond length 
is set to unity. 
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FIG. 1. (Color online) Real-space spin-spin correlation func- 
tion along the zigzag edge of ribbons. Points are the QMC 
and MFT data and lines are the fits with Eq. Inset: Ge- 
ometry of the graphene nanoribbon with zigzag edge, L = 6, 
ly = 4, and lattice size iV = 2xLxW = 48. The differently 
shaded circles represent the positive and negative values of 
the local magnetization (Sf) — ^{iii^-f ~ rii^i}, calculated from 
MFT ai U/t — 2. Radii are proportional to the magnitude. 



The Hamiltonian of the Hubbard model reads 



H 



(1) 



with 



nearest neighbors on the lattice, a —^,\r and 
rii = c[ci. We concentrate on half filling, i.e., the number 
of electrons is equal to the number of lattice sites. MFT 
for the 2D periodic honeycomb lattice yields a transition 
from a semimetal to an antiferromagnet at Uc/t = 2.23 
21] whereas quantum Monte Carlo (QMC) simulations 
reveal a metal- insulator transition for Uci/t = 3.5(1) and 
a further nonmagnetic insulator-antiferromagnetic insu- 
lator for Uc2/t = 4.3(1) [l^l. In this work, we focus on 
the interaction range U /t < 2 since we are interested in 
the behavior of the system in the weak-coupling regime. 
Estimates of the local Coulomb repulsion U for graphene 
12{ lie around U/t ~ 1 and thus fall indeed into this 
regime. 

The inset of Fig. [1] shows the MFT result for the local 
magnetization at the edge ofaVF = 4,i = 6 ribbon at 
U/t = 2 and illustrates the edge-state magnetism which 
is expected for weak Coulomb interactions U/t > [§- 
[i3|. This is reflected by the MFT result for the spin- 
spin correlation function which is shown by triangles for 
W = 6, L = 45, and U/t = 2 in the main panel of Fig. [T] 
the correlation function rapidly approaches a constant 
value as a function of distance x along the edge with 
very little dependence on the width W of the ribbon (not 
shown) . 

By contrast, the QMC results for the static spin-spin 
correlation function shown in the main panel of Fig. [1] 



exhibit a clear decay, at least for narrow ribbons with 
W = 2 (circles) and 4 (squares). Note that the rapid 
decay ior W — 2 is consistent with a previous density- 
matrix renormalization group study [17|. However, as 
the width increases, the decay becomes slower and the 
QMC result for 14^ = 6 in Fig. [1] (diamonds) is already 
qualitatively similar to the MFT result. 

For a more quantitative analysis we fit the QMC data 
for the longest available ribbons with a function combin- 
ing exponential and power-law behavior 



(SoS^) = C (a 



(^|ai| 



-(L|ai|-a;)/5 



(2) 

From a fit for L = 60 at 14^ = 2 we estimate ^/|ai| = 
4.0 ± 1 and = 0.5 ± 0.3 and at = 4 the L = 54 
data yield C/|ai| = 15 ± 4 and = 0.2 ± 0.15. At 1^ = 
6 it is no longer possible to distinguish the correlation 
length from infinity or the exponent 77 from zero for the 
largest available system size within our study {L = 45). 
Indeed, the W — 6 QMC data can equally well be fitted 
by predictions based on an infinite correlation length. We 
thus arrive at a picture similar to even-leg spin ladders 
[23! or integer spin- 5 quantum chains 24 [: ribbons with 
even W have a finite spin-spin correlation length 16[ , but 
the correlation length rapidly grows with W such that for 
practical purposes it can be considered infinite for W >6. 
In this sense edge-state magnetism is found already for 
ribbons of moderate widths. In order to check if MFT 
becomes also quantitatively accurate for wider ribbons, 
one would need to go beyond the system sizes which we 
can presently access by QMC calculations. 

Since the ferromagnetic behavior at the edges remains 
fluctuating in a proper treatment of the Hubbard model 
P^. we now turn to the dynamic properties. In this re- 
spect, the local single-particle spectral function A{uj) is 
a particularly important observable since it can be ob- 
served in scanning tunneling microscopy (STM) exper- 
iments [2^, [2^. However, to the best of our knowl- 
edge, previous theoretical studies of the interplay of 
edge-state magnetism and spectral functions of graphene 
ribbons are restricted to density-functional approaches 
111 . |27| with an inherent mean-field type approximation. 
In QMC we measure the momentum-resolved Green's 
function in imaginary time G{q,T), and then apply a 
stochastic analytical continuation scheme [2^ [2^ to ro- 
tate G{q,T) on the r axis to A{q,uj) on the w axis. Fi- 
nally, the local spectral function A{uj) is obtained by inte- 
gration over the momentum q along the ribbon direction. 

Figure [5] shows the single-particle spectral function 
both on the zigzag edge ^(a;)edge and inside the bulk 
^('^)buik for ribbon of width W = AatU/t = 1. The local 
spectral functions in the main panel were subjected to a 
Lorentzian broadening of Aw = 0.02 A(a;)odge exhibits 
a dominant low-energy peak at w/t « 0.1 which is absent 
in ^(a;)buik- This peak can be traced to a flat region of a 
single-particle dispersion visible in A(q, w)cdgG (compare 
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FIG. 2. (Color online) Spectral functions of a width W = 4 
ribbon at U/t = 1. QMC data are for a ribbon of length 
L = 30 whereas MFT corresponds to the thermodynamic 
limit. The inset shows QMC results for the momentum- 
resolved spectral function A{q,uj) along the zigzag edge, while 
the main panel shows the local spectral function A{uj) subject 
to a Lorentzian broadening Atj = 0.02 1. Note the prominent 
low-energy peak in A(aj)edge at u)/t ~ 0.1 which is absent in 

^('^)bulk. 



the inset in Fig. \^ . The agreement between MFT and 
QMC is excellent for this low-energy peak and the over- 
all agreement is also good at higher energies. The three 
additional sharp features in the QMC result for yl(a;)cdgc 
in the region 0.2 < uj/t < 0.6 can be traced to a finite- 
size momentum discretization effect and are reproduced 
by MFT if the latter is also restricted to a length L = 30 
ribbon. We conclude that the agreement between MFT 
and QMC for A{uj) is remarkable even at a quantitative 
level, at least for weak interactions U/t < 1. 

Having gained confidence in the accuracy of the MFT 
results for A{lo), we use it to analyze a bigger and 
more realistic system with W = 48, corresponding to 
a graphene nanoribbon of 10 nm width. The evolution of 
the low-energy peak on the edge with the Coulomb repul- 
sion is shown in Fig.[3^a). The energy Wmax correspond- 
ing to the maximum intensity of the spectral function 
increases with the Coulomb repulsion U and is located 
at Wniax = only for U = 0, i.e., the noninteracting sys- 
tem. In combination with the fact that this peak exists 
only on the ferromagnetic edge (compare Fig. [2]) we con- 
clude that it is a clear dynamic signature of edge-state 
magnetism. Figure [31[b) demonstrates that Wmax is in- 
sensitive to the actual width of the ribbon and confirms 
that MFT is accurate for U/t < 1 (the fact that MFT 
becomes less accurate as one approaches the mean-field 
phase transition at Uc/t = 2.23 [2l| is not surprising). 

The momentum-resolved spectral function [shown in 
Fig. |3Uc) for U/t = 2] reveals that the feature at a;,„ax 
can be traced to the maximum of a single-particle band 





0.1 0.2 0.3 0.4 0.5 4 - |^ 2ir 

uj/t q/\ai\ 

FIG. 3. (Color online) (a) Single-particle spectral function 
yl(cj)edge for several values of U/t and a ribbon of width W = 
48. (b) Comparison of the position i^max of the maximum of 
A{uj)cdge between MFT for W = 4, 48 and QMC for W = 4 
and L — 30. (c) A{q, cj)cdgc for a ribbon with W = 48 at 
U/t = 2. Unless stated otherwise, results are based on MFT 
for long ribbons and a Lorentzian broadening Aoj = 0.02 1. 

at 9/|ai| = TT where large matrix elements and a van 
Hove singularity combine to yield a maximum intensity 
of the local spectral function. The true single-particle 
gap Asp is located in the vicinity of (7/|ai| = 27r/3 and 
47r/3. For the = 48 ribbon shown in Fig. |3] it is 
only Asp/t = 0.037 for U/t = 2, resulting in a fill-in 
of spectral weight in the local density of states below 
the "pseudogap" Wmax [compare Fig. El^a)]. The value 
of the single-particle gap Agp depends strongly on the 
width of the ribbon; for the W — A ribbon shown in 
Fig. [2] it is quite close to the pseudogap Wmax- Note 
that the data in Fig. [S^a) are subjected to a broadening 
Aw = 0.02 1. With a higher energy resolution, one would 
observe further features, in particular another van Hove 
singularity at Agp, albeit with a much smaller weight 
than at Wmax- 

Another probe to exhibit the peculiar magnetic prop- 
erties of zigzag nanoribbons is provided by the spin spec- 
tral function S{q,Lo) on the edge. From the enhanced 
ferromagnetic correlations along the zigzag edge at finite 
values of J7, one might expect S{q, w) to exhibit collec- 
tive spin-wave (magnon) excitations characteristic of a 
ferromagnetic background, with a quadratic low-energy 
dispersion relation near q ~ Q. However, as pointed out 
in Refs. 0, ll^l, the magnetic excitations are affected by 
the antiferromagnetic coupling between the two edges of 
the nanoribbon that is mediated by the bulk conduction 
electrons. These antiferromagnetic correlations result in 
a linear contribution to the dispersion at small wave vec- 
tors (see also Ref. [l3[). Furthermore, for energies enter- 
ing the particle-hole continuum, the collective modes can 
decay into strongly damped Stoner modes [2G] . 
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which has observed low-energy spectral features by STM 
and compared those to MFT for isolated chiral graphene 
nanoribbons. 



FIG. 4. (Color online) QMC results for the spin spec- 
tral function S{q,u}) along the edge of a zigzag ribbon with 
L = 30, IV = 4 and U /t = 2. The inset shows the linewidth 
(FWHM) of the dominant mode in S[q, u) as a function of 
momentum q. 



Our QMC data in Fig. U] indeed exhibit such a sce- 
nario for a ribbon with W = A: an approximately lin- 
early dispersing sharp mode is observed up to a wave 
vector g/jail « O.Stt with a threshold energy wth/i — 
0.42 ± 0.04, compatible with twice the corresponding 
single-particle gap of Agp/i — 0.2 ± 0.02. The inset of 
Fig.|3]shows the corresponding evolution of the linewidth 
(FWHM) of this low-energy mode that drastically in- 
creases beyond (z/|ai| w O.Stt. 

In summary, we have investigated the static as well as 
dynamic properties of graphene nanoribbons with zigzag 
edges. Upon closer inspection, true ferromagnetic long- 
range order at the zigzag edge [l-[l3 of a finite-width 



ribbon turns out to be absent |14| - |17l |. yet there are 
strong ferromagnetic correlations which are essentially 
long range but for the narrowest ribbons. As a next step, 
we have identified a dominant low-energy peak in the 
local spectral function as a dynamic signature of edge- 
state magnetism, thus providing a theoretical guide for 
further STM and spin-resolved STM experiments. The 
agreement between QMC and MFT for the local spectral 
function is remarkably accurate for moderate Coulomb 
interactions, justifying a description of realistic geome- 
tries within a MFT framework. In particular, we have 
demonstrated that the position of the dominant low- 
energy peak at the zigzag edge of a nanoribbon is con- 
trolled mainly by the Coulomb interaction U such that 
STM experiments can be used to deduce the appropriate 
value for graphene. Lastly, we have presented QMC re- 
sults for the dynamic spin structure factor at the edge of 
a width W = 4 ribbon and demonstrated the presence of 
an approximately linearly dispersing low-energy collec- 
tive spin-wave excitation [7] that decays upon entering 
the particle- hole continuum (2C] |. 
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